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Understanding the nature of the exotic XYZ resonances is one of the 
open problems in hadronic spectroscopy. Despite the experimental ef¬ 
forts, the structure of these particles still lacks of an accepted theoretical 
framework. We propose to use the Zc'^ —)■ rjcP decays as a possible discrim¬ 
inant between two of the most popular models: the compact tetraquark 
and the loosely bound meson molecule. We show that the predictions 
obtained within the two pictures are signihcantly different and therefore 
the proposed decay channel might shed some light on the nature of these 
states. 
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Figure 1: Spectrum of the charmonium and charmonium-like states. Left panel: black 
and blue lines are the observed and expected charmonia respectively [1]. Red lines 
are the observed neutral exotic states. Right panel: observed charged exotic states. 
On the right of each panel we also report the open charm thresholds. 

1 Introduction 

In the past decade, several experiments observed a set of charmonium-like resonances 
sharing similar “exotic” properties — see [1] for a review. The spectrum of these 
particles, the so-called XYZ states, is reported in Fig. In particular, the obser¬ 
vation of charged resonances decaying into charmonia is a compelling evidence for 
the existence of four-quark states. The most popular models proposed to describe 
the internal structure of these particles are the compact diquark-antidiquark (or just 
tetraquark from now on) I21E], the loosely bound meson molecule HE], the so-called 
hadro-charmonium [B] and the gluonic hybrid [7j. None of these models has been 
generally accepted as the right one yet. 

Despite the lack of observation of many of the states predicted by the tetraquark 
model, the recent discovery of two baryon-like resonances with opposite parities de¬ 
caying into J/ipp ^ has given new interest to the diquark models [9]. 

In this work we summarize what developed in ra. focusing our attention on two 
charged states: the Zc(3900) and the Z'(4020). The first one has been observed Iw 
BES and Belle in the J [TT] final state and by BES in the [DD*)^ [12] channejj 
while the latter one has been observed by BES in hcir^ and {D*D*)^ [T3]. For both 
of them the most likely quantum numbers are = (1+) !■*■“, where the charge 

conjugation refers to the eigenvalue of the neutral isospin partner. 

Similarly to what happens to many other exotic resonances, the Zc and Z' lie very 
close in mass to the DD* and D*D* thresholds respectively (see again Fig. and 
therefore have been interpreted as loosely bound molecules m, despite being slightly 

*Here and in what follows charged conjugate modes are undestood, unless otherwise specified. 
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above threshold. A state with mass around 3882 MeV is however also 

predicted by the constituent diquark-antidiquark model [2], together with its radial 
excitation around 4470 MeV. These states have been identihed respectively with the 
Zc{3900) and with the recently discovered Z(4430) [15]. While in the hrst version 
of the tetraquark model (the so-called “type 1”) the Z'(4020) was not included, it is 
now nicely accomodated by a recent “type 11” paradigm [3]. 

Here we show how the Zc'^ —)■ rjcp decay channel can be used as a tool to differen¬ 
tiate between two of the possible internal structures of these charged resonances. For 
the Zc(3900) —)■ rjcP process, some considerations about the consequences of the heavy 
quark spin symmetry within the molecular and tetraquark pictures have already been 
explored in [T6] . 

2 Compact tetraquark 

In the constituent diquark-antidiquark model the Hamiltonian that describes the 
interaction between the components of the state is given by: 

= ( 1 ) 

i a,ij^j 

where ruj are the masses of the constituents, Kij are unknown couplings. Si are spin 
vectors and A“ are the Gell-Mann matrices. In the type 1 model [2], the couplings 
are extracted from the spectrum of ordinary mesons and baryons and the mass of the 
X(3872) is used as an input to determine the diquark mass, mjcg]. In the type 11 [3], 
instead, since the spin-spin interaction is a contact term, it is assumed that the only 
relevant interactions are those within the diquarks and hence all the couplings are set 
to zero except for Kcq = Kcq- Depending on which of the two ansatz is chosen, the 
physical states X(3872), Zc(3900) and Z'(4020) are identihed with different combi¬ 
nations of the eigenstates of the Hamiltonian in ([^. 

As far as the decays of tetraquarks are concerned one can instead resort to the 
well-known heavy quark spin symmetry im to write the amplitudes for the decay into 
charmonia as a Clebsch-Gordan spin factor times a transition matrix element [10] . 
This is valid up to corrections of order Agci^/'^c — 25%, where rric — 1.5 GeV is the 
constituent charm quark mass. 

For the processes of interest, the most general Lorentz-invariant matrix elements 
that behave properly under parity and charge conjugation are 

{J/^{p,p) 7r(g) |Z(A,P)) =gzpn^-V, {vdp) p (e, g) |Z (A, P)) =gzn,pX-e, (2a) 

{K {p, r,) n{g)\Z (A, ^’)) = ^ (2b) 
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where A, rj and e are polarization vectors, p, q and P are fonr-momenta and the gs 
are effective conplings with dimension of a mass. 

Since we have no information on the conplings, some kind of ansatz is reqnired 
here as well. To test the degree of model dependence of onr calcnlation we make two 
different assnmptions: (a) we neglect the spatial dependence of the wave fnnctions and 
hence assnme that the conplings are nniversal, the differences between the different 
matrix elements being only of kinematical natnre; (b) we nse a dynamical model 
recently developed in [TH]. In this pictnre the conplings for the decays into charmonia 
are proportional to the propabability density of the charmoninm itself compnted at 
a distance rz- The last one is defined as the classical tnrning point of a diqnark- 
antidiqnark pair moving away from each other and interacting via a spinless Cornell 
potential. 

In Tab. [^we report the predictions obtained within the tetraquark model. 



Kinematics only 

Dynamics included 


type I 

type II 

type I 

type II 

BP (Zc -)■ rjeP) 

X 102 

n 41 +0-96 

(2.3l?i) X 10^ 

-0.17 

BP{Ze^ J/'lpTl) 

BP{Z'^^ hen) 

(1.2«;*) X 102 

6.61S‘ 



Table 1: Predicted ratios of branching fractions for the Zc'^ states according to the 
main tetraqnark models. The hrst and second colnmns are compnted nnder the 
assnmptions (a) and (b) respectively, as explained in the text. Both type I and type II 
models give the same predictions for the BTZ (Z' —)■ pcp) IBP- (Z' —)■ hcVr), since both 
he and rje have spin Sec = 0. The errors are estimated via a toy MC simnlation. 


3 Loosely bound molecule 

As already mentioned, in the molecnlar pictnre the Z^^ is interpreted as a D^*^D* 
loosely bonnd state. The interaction between the exotic particles and the heavy 
and light mesons is commonly described by means of the so-called Non Relativistic 
Effective Field Theory (NREFT) [5]. This is a non-relativistic limit of the Heavy 
Qnark Effective Theory (HQET) together with the Chiral Effective Field Theory 
(xEFT). The complete Lagrangian of interest for onr stndy is fully reported in [10] . 
together with the choice of couplings for the interaction between the different fields. 
The term describing the interaction between the Zc^ and the charmed mesons is given 
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Figure 2: Possible one-loop diagrams for the Z,. (upper figures) and the Z' (lower 
hgures) decaying into r]cP- The charge conjugate diagrams are omitted. 


by 


d') 






w 


+ h.c., (3) 

and Hia are the HQET helds for the doubly heavy ZcS states and for 


where Z 

I 

the D mesons respectively. The z^''^ are, instead, unknown effective couplings. See 
again ng for details and definitions. In principle, such an effective theory is a valid 
description of the decays of the Zc'^ regardless of its internal structure since the form 
of the interaction is only dictated by symmetry considerations. The molecular nature 
of a state is imposed by forcing it to couple to its own constituents only. Therefore, 
the decays into final states different from the latter ones (charmonia in our case) can 
only happen via heavy meson loops. The most relevant one-loop diagrams for the 
zP —>■ rjcP process are reported in Fig. [2 

Moreover, since the molecular states are near threshold, the velocities of the 
mesons inside the loops are given by vx — \/\Mz — 2 Md\/Md and are typically 
small. This allows to perform a power counting procedure in order to estimate the 
relevance of a certain diagram [5]. Using this technique, one hnds that the omission 
of diagrams with more than one loop introduces a 15% relative error on each single 
amplitude. 

Given the previous set up, the predictions obtained within the meson molecule 
framework are: 


BTZjZ^ Pep) ^ . „+2.5^ X 1 n-2 . Bll{Z'^ Pep) 

Bn{Z^^ J/tl)Tl) V • - 1 . 7 ; > i37^(Z' ^ heTT) 


(1.01°:®) X 10-". (4) 
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Figure 3: Likelihood curves for BlZ{Zc —)■ ricp)/B7l{Zc —)■ J/tjjTr) (left) and B1Z{Z'^ —)■ 
rjc p)/B7l{Z'^ —)■ hcTi) (right). The red curve is the molecular prediction, whereas the 
black one gives the predictions for the dynamical type I tetraquark model. The green 
(yellow) bands represent the 68% (95%) conhdence region. 


As an additional result one can also assume that the total width of the Zc'^ is 
saturated by the D^*^D*, pcp, hc7i, J/V'tt and '0(2*5') vr hnal states and therefore ht 
the couplings to the constituents from the experimental data. This gives 

I 1 = (l.260K;i^) GeV-^/2 and | | = (O.580)J:?^) GeV^/^ (5) 


Once these couplings are given one can also make the following predictions for the 
comparison between the two charged resonances decaying into the same hnal states: 


Bn{Z^ ^ ^cTt) _ ^ ^^+0.21 Bn{Z, J/07r) _ ^ ^^+0 49 
i37^(Z' ’ f37^(Z' ^ J/0vr) “ ' 


( 6 ) 


4 Conclusions 

We can now properly compare the predictions obtained within the two models pre¬ 
sented — see Fig. As one can see, according to the dynamical type I tetraquark 
model the Z^ ^ Pcp decay should be enhanced with respect to the already observed 
Zc —>■ The opposite is expected in the meson molecule picture and the two 

predictions are separated by more than 2a (95% G.L.). A similar thing holds for the 
Z'c ^ Pc P decay with respect to the Z'^ ^ he ti one. In the last case, however, the 
predictions for the type I and type II models are the same and hence the result is 
more model independent. The values obtained under the assumption of no dynamics 
for the tetraquark turn out to be even more separated from the molecular ones. For 
the Zc in the type II paradigm, instead, the two models give predictions which are 
compatible within 2a. 
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Lastly, the results reported in Eq. ([^ show that in the molecular picture one 
expects l3TZ{Zc —)■ hcn)/BlZ{Z'^ —)■ hcvr) < 0.88 and BTZ{Zc —)■ J/'ip tt)/B 7l{Z'^ —)■ 
J/ipir) < 1.86 at 95% C.L.. This means that the two charged resonances should be 
seen in both the hcvr and J/tpn hnal states with comparable rates. While this seems 
to agree with the data in the hrst case, where a small hint of Zc is seen, it might be 
at odds with the experiments for the J/xpm channel, where no Z' has been observed 
so far. 

In conclusion, we showed how the analysis of the Zc'^ VcP decay can be used as 
a probe of the internal structure of these charged states and hence provide a tool to 
discriminate between two of the most accepted models for the exotic XYZ mesons. 
Experimental data on this channel could therefore shed some light on the now long¬ 
standing question about the nature of the Zc and Z' resonances. 
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